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In 1983 Haldane [1] argued that there is a fundamental difference between 
integer and half-integer spin chains (for a recent review, see [2]). Certain quantum 
antiferromagnets with integral spin - including the isotropic Heisenberg chain - should 
have the following properties: 
(i) the groundstate is unique, 

(ii) there is an energy gap between groundstate and excited states, 
(Hi) groundstate correlations have exponential decay. 

This will be called 'Haldane scenario' in the following and is in contrast to the 
behaviour of isotropic half-integral spin chains which are expected to have no gap 
and algebraic decay of correlations. This scenario has been verified in a number of 
experiments [3-6]. 

The first model for which all these properties could be proven rigorously was 
introduced by Affleck, Kennedy, Lieb, and Tasaki [7]. This so-called VBS-model 
(Valence Bond Solid) and some generalizations of it have subsequently been studied 
intensively, see e.g. [8-13] and references therein. 

In this letter we present a large class of antiferromagnetic spin-1 chains for which 
the groundstate can be given explicitly as a product of matrices. We like to point out 
that our class differs from other models with product groundstates as the well-known 
Majumdar-Gosh type models [14] (for further references, see e.g. [7,12]). 

The most general spin-1 chain with nearest-neighbour interactions and the 
following symmetries 

(a) rotational invariance in the x — y— plane, 

(b) invariance under S z — > — S z , 

(c) translation and parity invariance, i.e. invariance under the exchange j <-> j + 1, 
can be written in the following form: 



L 




3 = 1 

h jtj+1 = a A 2 + a x (A 3 B 3 + B 3 A 3 ) + a 2 B 2 J + a 3 A 3 + a 4 B 3 (l + B 3 ) 



(1) 



+ a 5 ((S*) 2 +(S* +1 ) 2 )+c, 



with real parameters a 3 , a constant c, and periodic boundary conditions. The nearest- 
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neighbour interactions are 




(transverse) , 
(longitudinal). 



= ±(S?±iS*). Spin 



chains of the type (1) have been studied extensively in recent years using various 
analytic and numerical methods [15-22]. Only a few special cases of (1) are known to 
be exactly solvable (for some references, see [7,12,13,19]). 

In the following the constant c in (1) will be adjusted to have the groundstate 
eigenvalue of /ijj+i at 0. Thus we have 



i.e. all eigenvalues of 7i are non- negative. 

We shall be interested in the antiferromagnetic case of the model where the 
groundstate is characterized by S% otal = 0. In the following we show that in a 
certain subspace of the -parameter space the antiferromagnetic groundstate can 
be constructed in the form of a matrix product of single-site states [13]. These states 
will be called Matrix-Product-Groundstates (MPG) in the following. In [13] we have 
shown that the groundstate of the VBS-model and its (/-deformed generalization [19] 
is a special MPG. 



Using the S z eigenstates \0)j and |±)j we define a local 2x2 matrix at each site 

3 b y 



with nonvanishing parameters a, a ^ and the global (antiferromagnetic) state 



Here '<g>' denotes usual matrix multiplication of 2x2 matrices with a tensor product 
of the matrix-elements. We shall show in an extended paper [23] that the ansatz (4) is 
most general for the antiferromagnetic region of (1) in which we are mainly interested 
here. That |V>o) is antiferromagnetic, i.e. S£ otal \ipo) = 0, can be checked easily. 

We demand that the MPG \ipo(a, a)) is groundstate of H with eigenvalue 0, which 
requires certain conditions to be satisfied. Since we have the obvious implication 




(3) 





\^ (a,a)) = Trace g 1 <g> g 2 <8> ... <S> 9l ■ 



(5) 




(6) 
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and we have the product property (5) it is sufficient to show that 



hj,j+i(gj®gj+i) = 0, (7) 

i.e. each of the four elements of gj(B)gj + i which are linearly independent is a groundstate 
of hjj+i. Calculating the product in (7) explicitly one can show that (3), (7) are 
satisfied provided the following equalities 

1) a = signa 3 , 2) aao = as — a±, 

(8) 

3) a$ = I CK3 1 + «o(l — a 2 ), 4) a 2 = olqcl 2 — 2|a 3 |, 

and inequalities 

a ^ , a 3 ^ 0, 

(9) 

ct4 > 0, ao > 

hold. The inequalities guarantee that all other eigenvalues of fajj+i are strictly 
positive. Thus \ipo(a,a)) is groundstate of the model (1) with eigenvalue if only 
the restrictions (8) are satisfied. In addition, if also the inequalities (9) are satisfied 
this groundstate is unique for any chain length L and there exists a finite gap A > 
to the excitations (in the thermodynamic limit L — > oo) as will be shown in detail 
in [23]. The first assertion can be understood from the fact that the MPG (5) is 
an optimal state in the sense that it is the product of all the local groundstates of 
hjj+i. The second assertion can be proved by a slight generalization of the proof [7] 
for the VBS-model which is recovered as the special case a = 2, a = 1, as = 3«o > 0, 
ct2 = — 2cto and a^ = 3ao [7,13]. Thus (i) and (ii) of the Haldane scenario are verified. 

It is clear by continuity that on the boundaries (equality signs in the inequalities 
(9)) \ipo(a,a)) is still groundstate, but no longer unique. 

In order to verify point (Hi) of the Haldane scenario we calculate the groundstate 
correlations explicitly. This can be done most easily using the transfer-matrix method 
of [13] and is explained in detail in [23]. We find for the longitudinal 2-site correlation 
for L — > oo and r > 2 

(s^) = --4-^(^) r (io) 



(1 - |a|) 2 \1 + \a 

and for the transversal correlation 



Ws-;:)= -M [a + sigH ( ) . (in 



Furthermore we find 



( 5 ;**> = o, ((^) 2 ) = -^. (12) 

Obviously the correlations (10), (11) decay exponentially with longitudinal and 
transverse correlation lengths 

C^Ml + lal). (13) 

The string-order-parameter [24-26] a z = lim^oo ^— Sf exp (itt X^=2 ^l) &r ^ can 
also be calculated easily. It is simply given by ((S^) 2 ) and thus a z ± 0. 

Now we briefly discuss the structure of the solution manifold S defined by (8) 
and (9) but refer to [23] for a detailed discussion. Apart from the additive constant 
c the model (1) is defined by six parameters ao, ■ ■ ■ , a$ including a trivial scale. Two 
of the four conditions (8) are satisfied by fixing the two free parameters a, a of the 
MPG (5), the remaining two conditions, i.e. 3), 4) of (8), reduce the number of free 
parameters of the model to four. Thus S is a 4-dimensional parameter subspace. The 
parameters may be chosen as ao, a^ > 0, as ^ and a\ such that a ^ 0. Obviously 
S has external boundaries, ao = and a^ = 0, and two internal boundaries, namely 
«3 = and a = 0. 

At «4 = the MPG is still a groundstate but it becomes degenerate with the 
ferromagnetic eigenstates ! + •••+) and | — •••—) which both have gapless spin- wave 
states adjacent to them (see also [23]). We have a first-order phase transition to the 
ferromagnetic regime at a^ = 0. 

At ao = the parameter a is not uniquely determined by (8) but \ifJo( a , c)) spans 
a L/2-dimensional space of degenerate states. One of these states is |0 - - - 0) which 
again has gapless spin-wave states adjacent. Here we have a first-order transition to 
an (antiferromagnetic) regime. 

Within S we have two internal boundaries as = and a = where also phase 
transitions occur. At as = a = sign(«3) jumps from + to — indicating a first-order 
transition. This is seen explicitly in the behaviour of the transverse correlation (SfS^) 
(11). The two MPGs \ipo(a, +)) and ^(a, — )) are degenerate with each other and also 
with the states | ±0±0 • • •) and |0±0±- ■ ■), indicating a vast groundstate degeneracy. 
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er 1 = ^ 



i -+- a 



1 - a 



At the other internal boundary a = the MPG becomes \ipo(Q,a)) = |0 • • - 0) , 
i.e. all spins are in the x — y-plane. As in the case cto = we find low-lying spin- 
wave excitations and a vanishing gap. As is obvious from (13), the correlation lengths 
diverge with a — > 0, thus we have a critical transition. 

In conclusion we note that we have found an interesting realistic model for which 
the groundstate can be determined exactly and which shows a number of different 
phases. 

We also like to mention that a generalization of the above considerations to 
arbitrary spin S - including non-integral S - is possible [23]. Also, using MPGs as 
variational states for the bilinear-biquadratic spin-1 chain (i.e. the isotropic case of 
(1)) gives excellent results for the groundstate energy [23]. 
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